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The FLUKA calcnlation of the atmospheric nentrino flnxes have been cross-checked by 
comparing predictions on lepton flnxes in atmosphere to experimental data. The depen¬ 
dence of predicted nentrino flnxes on the shape and normalization of primary spectrnm 
is also investigated. 


1. Introduction 


The two essential featnres of the atmospheric nentrino fluxes, namely the up-down 
symmetry and the /i/e ratio, are now considered to be largely independent from the 
calculation scheme[^. Theoretical uncertainties should still have an impact, together 
with experimental systematics, in the determination of the relevant physical parameters 
of neutrino oscillations[^. However, the Super-Kamiokande allowed region in the Am? — 
sm^(20) space remains more or less invariant [Q when the Monte Garlo predictions change 
from the original HKKM resultsto the FLUKA onesf^, despite the fact that the two set 
of predictions have important differences. The most signihcant one is the normalization 
of the flux, which, for the same neutrino cross-sections, gives a 15-^20% lower event rate 
in the case of FLUKA. The same ratio approximately exists between FLUKA and the 
original Bartol flux|^ for the same primary flux. This difference is due to the features 
of the hadronic particle production model. Our claim is that FLUKA has already given 
a more accurate prediction than the previous calculations. Also in the recent work of 
the Bartol group |^, in which some feature of their hadronic interaction model have been 
corrected, a neutrino flux very close to our predictions is now obtained for high cut¬ 
off sites, like Super-Kamiokande. It is likely that in the case of Super-Kamiokande 
analysis the differences between fluxes are in practice obscured by the assumed systematic 
uncertainties. However, in view of possible future improvements, it is important to clarify 
the differences among the existing predictions. The effects due to the 3-D treatment of 
the shower, hrst presented in ref. |§], have been reproduced by other groups P,pR| and are 
by now accepted as the correct prediction. The uncertainties in the neutrino rates are 
instead less clear, but they can be attributed to three independent sources: i) the primary 
cosmic ray spectrum, ii) the hadronic interaction model and hi) the neutrino-Nucleus cross 
sections. The hrst goal of this paper (section is to provide further evidence that the 
particle yields predicted by the FLUKA model |P,|l| are in the right range. As a second 
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step, we consider the dependence of the predicted v flnx on the choice of primary spectrnm 
(section |). The fnndamental topic of nentrino cross-sections is instead postponed to a 
fntnre work. In the conclusions we summarize the arguments, advancing some criticism 
about other new attempts of v flux calculations. 


2. New benchmarks for the FLUKA predictions 


FLUKA is constructed upon a set of theoretically inspired models of hadronic interac¬ 
tions, adjusting parameters only on the basis of results coming from accelerators. Cosmic 
ray data can however be used to check their quality. Recently two remarkable results 
have been achieved in this held. The hrst one is the reproduction of the features of the 
primary proton hux as a function of geomagnetic latitude as measured by AMS|]^, thus 
showing that both the production of secondary nucleons and the geomagnetic effects and 
the overall geometrical description of the 3-D setup are well under control. The same 
work also shows that the huxes of secondary and e~ measured at high altitude are 
reproduced, and this is instead more directly linked to the yield of mesons produced in 
primary interactions. The second achievement is the good reproduction of the muon data 
in atmosphere as measured by the CAPRICE experiment ||12||, both at ground level and 
at different floating altitudes [|T^, when starting from the same primary flux (Bartol ht)P| 
used to generate atmospheric neutrinos. The fluxes of atmospheric muons are strictly re¬ 
lated to the neutrino ones, because almost all z/’s are produced either in association, with, 
or in the decay of /i^. The level of agreement reported in ref.]^, without any adjust¬ 
ment of the model, allows to conclude that, the normalization of the v fluxes predicted by 
FLUKA is in the right range. This is true for a choice of primary spectrum constrained 
by the more recent data, in particular from AMS[|T^] and BESSfl^. The quoted primary 
flux of ref. IP is already in good, although not yet optimal, agreement with these data. 
Furthermore, the agreement exhibited by the FLUKA simulation separately for muons 
of both charges gives conhdence on the simulated yields of 7r+ and 7r“ production in air. 
This is one of the main sources of uncertainty in the prediction of the /x/e ratio. The 
decay chain of eventually results in Ue, while the decay of 7r~ generates Vf,. Due to the 
different interaction cross section this affects the experimental event rates. 


3. Neutrino Fluxes and Primary Spectrum 


For a given shower model, the level of agreement in a comparison to measured particle 
fluxes in atmosphere depends on the parameterization of the primary spectrum. The 
primary choice originally used in ref.p would force in FLUKA a too high normalization 
of the flux, since it is includes the data points taken from the compilation of Webber 
et al.|T^, now excluded by recent experiments. On the contrary, a constraint based on 


the primary data from CAPRICEwould not allow in FLUKA a satisfactory level of 
agreement with muon fluxes. As already mentioned, the Bartol £t used so far is not com¬ 
pletely satisfactory. The scientihc community is trying to establish a common reference 
for the primary flux to be used in all computations. An attempt in this direction has 
been proposed in ref.[|T^], where the statistical errors of measurements and the systematic 
differences between experiments are used to establish a band of uncertainty around the 
average result of this new £t. Contributions heavier than primary protons are also con- 
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sidered. A similar attempt has been proposed in ref[0. We have used the new Bartol 


£t of primaries to recalculate the FLUKA u fluxes. As an example, in Fig. we plot the 
ratio of the computed flux for the site of Super-Kamiokande as resulting from the new 
Bartol £t with respect to the old one, as a function of v energy. 
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Figure 1. Ratio of the computed flux at Super-Kamiokande as resulting from the new 
Bartol £t with respect to the old £t, as a function of u energy. The uncertainty band is 
also shown. 


The new parameterization brings to an increase of neutrino fluxes in the Sub-GeV 
range, while a depletion occurs in the Multi-GeV range. This behavior is still maintained 
when event rates are considered. Integrating in energy, the average prediction for the 
atmospheric neutrino event rates increases, approximately for all flavors, by ~5% with 
respect to the original evaluation, with an uncertainty (only due to the primary spectrum) 
of ~ i4%. 
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4. Conclusions 

The absolute normalization of atmospheric neutrino fluxes is still a subject of investiga¬ 
tion. The results summarized in this paper support the conclusion that the normalization 
of the FLUKA results should be close to reality and that the difference with respect to 
the original Bartol and HKKM results cannot be attributed to an insufficient particle 
production yield in the shower development in the FLUKA model. The knowledge of 
the primary spectrum is still an important issue, but the convergence towards a common 
reference, constrained by the AMS and BESS data allows to reduce the related uncer¬ 
tainty (± ~5%). After this study we can also conclude that the agreement of FLUKA 
predictions to Super-Kamiokande|^ data will certainly improve if the updated param¬ 
eterization of the proton and He fluxes are used. This should be more evident for the 
Sub-GeV events. This does not mean that the questions about the absolute values of 
neutrino cross-sections and the hadronic interaction descriptions are answered. In par¬ 
ticular we remind that hadron-Nucleus and Nucleus-Nucleus interactions in the energy 
range 1-^30 GeV are still insufficiently known to allow the construction of a fully reli¬ 
able model. In any case it seems difficult that the overall uncertainties are larger than 
15%^20%. Therefore, in our opinion, the results on neutrino flux of ref.[^, starting from 
a primary spectrum very similar to the one here considered, seem too low in the Sub-GeV 
region, not being compatible with the corresponding muon fluxes in atmosphere, unless 
significant systematic error in the data at high altitude are assumed. Furthermore, after 
our tests on direct Nucleus-Nucleus interaction, by means of an interface to the DPMJET 
code [pm 


we cannot share their conclusion about the existence of a significant bias intro¬ 


duced by the use of the superposition model. As far as the v fluxes proposed in ref.[pl| 
we believe that they are probably affected by some severe error (see also ref. [Q]). Beyond 
other considerations, those results would imply, without any supporting evidence, that 
either the experimental data on atmospheric neutrinos are strongly biased or that the 
present knowledge of neutrino cross sections are wrong by an unexplained large factor. 
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